arXiv: 1505.0001 lvl [hep-ph] 30 Apr 2015 


Accelerating the Discovery of Light Dark Matter 


Eder Izaguirre, 1 Gordan Krnjaic, 1 - 0 Philip Schuster, '[j] and Natalia loro 1 [-J] 

] Perimeter Institute for Theoretical Physics, Waterloo, ON N2L 2 Y5, Canada 
(Dated: May 4, 2015) 

We analyze the present status of sub-GeV thermal dark matter annihilating through Standard Model mixing 
and identify a small set of future experiments that can decisively test these scenarios. 


INTRODUCTION 


GEV-SCALE DARK MATTER 


That dark matter (DM) is a thermal relic from the hot 
early Universe is an inspiring possibility that motivates non- 
gravitational interactions between dark and ordinary matter. 
The canonical example involves a heavy particle interacting 
through the weak force (WIMPs). This scenario has moti¬ 
vated searches for DM scattering in underground detectors, 
for DM annihilation in the cosmos, and for DM production in 
high-energy colliders. These efforts achieve broad and pow¬ 
erful sensitivity to DM with mass between a few GeV and the 
TeV scale. 

A thermal origin is equally compelling — and, in simple 
models, predictive — even if DM is not a WIMP. DM with 
any mass from an MeV to tens of TeV can achieve the correct 
relic abundance by annihilating directly into Standard Model 
(SM) matter. However, the lower half of this mass range can¬ 
not be fully explored using existing strategies - an unfortu¬ 
nate situation that jeopardizes the legacy of the DM search 
effort. In particular, DM-nuclear and DM-electron scattering 
searches lose sensitivity precipitously for DM lighter than a 
few GeV or DM that scatters inelastically; limits on DM anni¬ 
hilation at low temperatures (most notably from the CMB) are 
irrelevant to many scenarios; and missing energy searches at 
high-energy colliders are blind to the interactions responsible 
for the DM abundance. 

In this paper, we sharply define the challenge of testing sub- 
GeV thermal DM. Below a few GeV, avoiding DM overpro¬ 
duction requires that a comparably light particle mediate DM 
annihilation mm. We focus on the case where DM annihi¬ 
lates through an s-channel mediator directly into SM state 
We classify annihilation mechanisms arising from dark-sector 
mixing with SM fields, define minimal models representative 
of each mechanism, and compute all relevant constraints on 
each model. We introduce a framework to compare all such 
constraints to one another and to the milestone provided by 
thermal DM freeze-out. Finally, we identify a small set of 
“flagship” experiments with complementary sensitivity — us¬ 
ing direct detection, /j-factory mono-photon, and fixed-target 
missing momentum strategies — that together can decisively 
test light DM annihilating through SM mixing. 


1 Sufficiently heavy DM can also annihilate into two mediators — a far less 

predictive scenario. 


Viable light thermal DM scenarios can be classified by the 
spins and masses of the DM and mediator, by whether the 
thermal abundance or a primordial asymmetry dominates the 
DM density, and by the mediator’s interactions with both DM 
and SM matter. SM symmetries substantially restrict the lat¬ 
ter interactions: vector mediators can mix with the photon or 
weakly gauge a SM global symmetry, while scalars can mix 
with the Higgs (or have axion-like couplings in extensions of 
the SM). Rare /i-meson decays largely exclude the scalar me¬ 
diator scenarios for sub-GeV thermal DM (see below), so we 
focus here on the possibilities for DM coupled through a vec¬ 
tor mediator. 

One concrete example is a scalar QED model of DM, where 
the “dark photon” A' is massive with coupling go = 1 /47t ap 
to scalar DM currents + c.c.. The DM scalar 

if and dark photon A' have masses m v and m A r respec¬ 
tively. The leading SM coupling to this dark sector allowed by 
symmetries is photon-A' kinetic mixing, C m i X = eF'^F^ 
where F' and F are the A' and photon field strengths, re¬ 
spectively ®0. After diagonalizing away the kinetic mixing 
term, the low-energy Lagrangian that describes dark-visible 
interactions is 


tint — A!^(te3^ M +gcJn), (1) 

where is the SM electromagnetic current. 

For this representative case, we can now ask what parame¬ 
ter ranges achieve the correct tp relic density. For niA' > m^, 
the rate of annihilations ptp* —> ff determines the relic den¬ 
sity. Neglecting m //m v corrections, the tree-level annihila¬ 
tion cross section at relative velocity /; re i <C c is 

_ 8 tt e 2 aa D mlv^ 

CTWrcl 3 (m\, — ) 2 + m 2 A ,T 2 ’ 


where T is the A' width. In the limit rriA> 3 > rn v , T, this 
cross-section depends on dark-sector parameters only through 
the DM mass m v and the dimensionless combination 


y = e 2 a D 


( rn-tp 

\m A ’ 
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(3) 


so matching the ip relic abundance to the observed DM density 
essentially fixes y as a function of m v (models with larger y 
can give rise to a subdominant component of the DM). Of 
course, near the fine-tuned region m A ' ~ 2m v , the precise 
milestone differs from that inferred from Eq. Q. 
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Pseudo-Dirac Thermal Relic DM 



Asymmetric Dirac Fermion DM 


Scalar Thermal Relic DM 




FIG. 1. Constraints and projections for representative vector-portal DM scenarios. For definiteness, we evaluate all constraints for 
uidm/wa' = 1/3 and (except for the LSNDxSIDM bound - see below), ctn = 0.5, near the perturbativity limit. The relic density, 
CMB, and direct detection contours scale roughly as e 2 a£>(mDM/mA') 4 (plotted on the t/-axis), and so are insensitive to separate factors in 
the above. For other constraints, this choice is conservative, in that smaller choices of and/or moM/m-A' would shift the shaded regions 
downward (see text); arrows denote the shift in sensitivity for old —> 0.05. We illustrate these constraints for (left) pseudo-Dirac/inelastic 
fermion thermal-relic DM, with splitting 5 > 100 keV, (center) asymmetric Dirac fermion DM, and (right) scalar elastic-scattering thermal 
relic DM. Dirac fermion thermal-relic DM is fully excluded by the CMB constraint and inelastic or asymmetric scalar DM is quite similar to 
the right figure, but with CMB and direct detection constraints weakened. CMB, self-interaction (SIDM), and direct detection constraints all 
depend on the xiv) abundance, and are computed assuming the full DM abundance, not the thermal abundance expected for given masses and 
couplings. In all plots, gray shaded regions (color online) represent traditional DM constraints (e.g. direct detection), while non-traditional 
accelerator probes are shaded beige. We note that pseudo-Dirac limits are modified (and new dedicated searches are possible GDI) if 8 is large 
enough that \+ can decay on detector length-scales. 


Before comparing existing data to this milestone, we com¬ 
ment on obvious and important variants of the model above. 
First, the DM may be a fermion instead of a scalar. A Dirac 
fermion x = (xi ■ X 2 ) (decomposed here into Weyl spinors) 
can couple to the A! through vector and/or axial currents. 
The axial piece leads to velocity-suppressed p-wave annihi¬ 
lation with scaling similar to Eq. while the vector current 
= Xi<f M Xi — X 2 ^ m X 2 leads to s-wave annihilation, and 
typically dominates. For this reason, we shall focus on the 
pure vector coupling. 

If the global symmetry under which xi ,2 have opposite 
charges is broken (e.g. by a higgs field that gives mass to the 
A'), operators such as Cbreak = <5xiXi yield mass eigenstates 
X± = l/v / 2(xi i X 2 ) s pht in mass by 5, with off-diagonal 
A' couplings Ci n t = A^x+d^X-- This exemplifies the in¬ 
elastic or pseudo-dirac scenario CD. Analogously inelastic 
interactions can also arise in the scalar case. 

Finally, for either scalar or fermionic DM, its total abun¬ 
dance may be set by a primordial particle/anti-particle asym¬ 
metry that dominates over the thermal relic abundance. In this 
case Eq. Q sets a lower bound on the collective interaction 
strength so that the symmetric component is sub-dominant. 

Each scenario above has a counterpart where the A' couples 
to a global symmetry current of the SM (e.g. baryon minus 
lepton number), rather than via kinetic mixing. The results 
that follow rely mainly on the A' coupling to electrons, and so 
apply equally well (with 0(1) corrections to the thermal relic 
curve) to these scenarios, unless the A! gauges a symmetry 
under which electrons are neutral, such as p — t number EH 

23 - 

Scalar Mediators To illustrate the strong meson-decay 


constraints on scalar mediators, we consider one explicit 
model: a scalar mediator $ that mixes with the Higgs boson 
and couples to a DM fermion x, with £ int = JT + 

y y <\>XX, where y x and are free parameters and the y, 
are SM Higgs Yukawa couplings, y t = y/2rrii/v with v = 
246 GeV. Such a $ can mediate the partly invisible //-meson 
decays B+ -A AT+(*)$ -A AT+W XX, with a rate computed 
(for on-shell <f>) in 03 23 - When to$ > M b — M K 
rn x , this process (with off-shell <f>) has similar kinematics 
to B + -A A+W vu, the limit on the latter CD implies 
y\yje%/m% < 1.6 x 10 -6 GeV -4 . The DM annihilation 
rate scales similarly, but with y t replaced by the much smaller 
electron and muon Yukawas. This bound rules out thermal- 
relic DM for m x < GeV. The limits for lighter <1» and on 
scalar DM are even stronger, and constraints on axion-like 
couplings to Standard Model matter are comparable within 
0(1) factors. We defer a complete discussion of these scenar¬ 
ios for future work Qo). 


EXISTING DATA CONFRONTS LIGHT DM 

Returning to the representative scenarios with a vector me¬ 
diator, we now assess how well they are constrained by cur¬ 
rent data. Fig. [T] quantifies each constraint in the plane of 
y = e 2 an (m^/mA 1 ) 4 vs. m,, (or similarly for a fermion x), 
to facilitate comparisons with the relic abundance target. The 
scalings below apply for toa' > 2m v ( x ), where the A' de¬ 
cays invisibly into tp (x) pairs, but the same experiments also 
constrain tp (x) production through a lighter off-shell A!. 
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Pseudo-Dirac Fermion DM (Visibly Decaying A ') 



FIG. 2. Same as Fig[T] but assuming = l.lm x , where the 
relic abundance is still achieved through s —channel annihilation, but 
A' decays visibly. Dashed blue contours show projected sensitivities 
from (291 for future A 1 —» t + l~ searches in the next few years. De¬ 
pending on the x and A' masses, models below the relic density line 
may still achieve viable thermal abundances through other processes 
— see (30) . 

CMB and BBN Although DM annihilations freeze out be¬ 
fore the era of recombination, residual annihilations can re¬ 
ionize hydrogen and distort the high-f CMB power spectrum 
HMD- These data constrain the total power injected by DM 
annihilations ED, which scales as the DM annihilation cross- 
section (and hence approximately as y) times a correction fac¬ 
tor calculated in 820) . The resulting constraint on y defini¬ 
tively rules out thermal-relic Dirac fermion DM (not shown in 
Fig. []}, but not the other scenarios. For scalar DM (Fig. [I] 
right), p -wave suppression of annihilations at late times re¬ 
sults in a weak upper bound on y\ the bound is likewise expo¬ 
nentially weakened for inelastic DM as the excited state can 
decay or be thermally depopulated before recombination G2 
(not shown in Fig. [TJleft). For asymmetric fermion DM (Fig. [I] 
center), CMB constraints imply a minimum value of y needed 
to depopulate the symmetric DM component m. Weaker 
bounds arise from the 3 He abundance during BBN 8231426) . 

Direct Detection Light DM scattering at halo velocities 
induces nuclear recoils below current experiments’ energy 
thresholds, but electron-scattering of sub-GeV DM (whose 
rate scales with y) is constrained by XENONIO data I27li281. 
For pseudo-Dirac DM (Fig. [I] left), tree-level scattering is in¬ 
elastic and kinematically forbidden for 8 > keV; elastic scat¬ 
tering arises from a one-loop box diagram scaling as y 2 . 

B-Factories Following [37] |32) . we reinterpret the BaBar 
search for an (untagged) T(35') —► 7 + invisible 133) to con¬ 
strain the continuum process e + e _ —► 7 + A'^ —> "/XX- 
For on-shell A' , the event rate scales as e 2 , independent of 
ct£>, and insensitive to niA> and m v in our kinematic regime 
(except for a shelf at toa' ~ 1 GeV due to a peaking 
background). To compare to the relic density target in the 


y = e 2 ctD vs. m v plane, we must fix these 

other parameters. To do so conservatively, we use the bench¬ 
marks rn^/ruA’ = 1/3 and old = 1/2 (smaller ratios and 
smaller ajj would overstate the exclusion of thermal relic DM 
by these results, while larger ratios qualitatively change the 
physics and larger a /> can run towards the non-perturbative 
regime 1341 ). It can be argued that ~ tua 1 is most nat¬ 
ural in any case, and the resulting constraint only improves 
linearly with decreasing old on the y vs. m v plane. 

LEP and LHC Electroweak precision tests at LEP con¬ 
strain the existence of a new massive photon. In particular, ki¬ 
netic mixing induces a shift in the mass of the Z° boson, and 
the constraint depends on e and only mildly on toa' ■ In the 
limit rriA' ro^o, e is constrained to be e < 0.03 1351136 1. 

DM could also give rise to missing energy Q£t) signals at 
the LHC, for example in association with a jet or photon. We 
recast the results of a CMS DM search | [37) in the monojet 
-\-$t channel using 20 fb -1 at y/s = 8 TeV, which bounds 
e < 0.1. The signal yield for this search scales as e 2 and 
is set by kinematic cuts, whose efficiency is roughly flat for 
the MeV-GeV masses, so the //-scaling is similar to that of 
/j-factory bounds. 

Beam-Dump Experiments Fixed-target beam dump ex¬ 
periments offer another powerful probe of light DM. LSND 
data sets the strongest constraints for mA> < m n 0 GHffiD, 
via pN —> 7 T° + X, with 7 r° —► 7 A' —> 7 xx, with subsequent 
DM scattering off an electron, nucleon, or a nucleus mimick¬ 
ing neutrino neutral-current scattering. The electron beam- 
dump experiment El37 is sensitive to eN —> eNA', A' —» 
XX with subsequent DM-electron scattering in a downstream 
detector 8421 . In both cases, the production yield scales with 
couplings as e 2 and the detection probability as e 2 ao, so that 
overall yield scales as e 4 a^. The benchmark old = 1/2 
yields a conservative bound on y that improves as 1 /y/cw 
for smaller old- The itla' -dependence of these experiments’ 
yields is more complex, but again the choice m v /mA> = 1/3 
is conservative, with stronger constraints for larger ratios. 

DM Self-Interactions Constraints from the Bullet Cluster 
and from cluster lensing j43l [44) bound the self-interaction 
cross section to be cr S eif/rtiDM < cm 2 /g, which implies 
old < 0.06(MeV/m x ) 1//2 (mA'/10MeV) 2 and is only rele¬ 
vant at low mass. Saturating this limit on old instead of using 
our benchmark old = 1/2 leads to a tighter constraint on y 
from LSND for low (m x < 10 MeV) masses, which is plotted 
as a green dashed line in Fig. |T|(again assuming a conservative 
benchmark m^/mA' = 1/3). 

Supernovae The production of any free-streaming xx pairs 
inside supernovae (SN) is constrained by observations of SN 
1987A 05] 06). As discussed in El, the SN core luminosity 
can be appreciable, but the scattering of x off baryons with 

9 T 2 

cross section cr = Attololdo - s £ l , where T$n ~ 10 MeV 

m A' 

is the core temperature, can be large enough for the \ diffu¬ 
sion escape time to fall below a few seconds. This occurs for 
a_oe 2 > O(few) x 10” 14 (or y > 10 16 in the above Figures) 
in the mA' ~ m x ~ 10 MeV range of relevance for these 
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reactions. Thus, SN constraints do not probe the thermal dark 
matter target near y > 10~ 15 at low masses, though a more 
detailed analysis of this physics is called for. 

Visibly Decaying Mediators If mA' < 2m v ( x ), the A' 
signals at colliders are quite rich — the accelerator-based con¬ 
straints discussed above have counterparts proceeding through 
an off-shell A! , and searches for visible decays A! —> £ + £~ 
become quite powerful. Numerous experiments 0 SEMI 
constrain this parameter space, and several more are expected 
to run in the next few years ll62l - i67l . These results and new 
experiments are summarized in li29l . The interplay of these 
constraints is illustrated in Figure [2] for inelastic fermion DM 
with representative parameter choices ijia' = lAm x and 
«£> = 0.5. The situation is qualitatively similar for the scalar 
and asymmetric DM scenarios. 

THE STATUS AND FUTURE OF LIGHT DM SCIENCE 

Present Status While scalar mediated light DM scatter¬ 
ing is tightly constrained by meson decay, vector mediated 
scattering is viable over a wide range of mass and couplings. 
Figs. □ and [2] illustrate that current experiments are several 
orders of magnitude away from decisively testing thermal 
vector-portal dark matter. Electron-scattering limits are 2- 
3 orders of magnitude short of the thermal target even at 
the optimal DM masses for these experiments; CMB energy 
injection bounds exclude the thermal target for symmetric 
Dirac fermion DM, but are far from constraining pseudo- 
Dirac, asymmetric, and scalar DM scenarios. Low-energy 
accelerator-based experiments have made significant progress 
exploring this parameter space - BaBar data does robustly ex¬ 
clude scalar vector-portal DM with few-GeV mass - but they 
too leave 2-A orders of magnitude of unexplored parameter 
space below a GeV. 

Future Probes With Existing Strategies Progress search¬ 
ing for light DM based on existing experimental strategies will 
be driven by three fronts: direct detection, B-factories, and 
fixed-target experiments. Direct detection through electron 
scattering is currently background-limited; future improve¬ 
ments depend on this relatively new field’s success in push¬ 
ing experimental thresholds to lower energies and minimizing 
backgrounds [ 27,, 281. In parallel, nuclear-recoil direct de¬ 
tection experiments are expected to lower their target masses 
and recoil energy thresholds enough to start exploring the sub- 
GeV region 116814701 . Belle II, an upcoming high-luminosity 
B-factory in Japan, can also significantly improve on current 
B-factory sensitivity to light DM if it is instrumented with a 
mono-photon trigger l32l . Together, these experiments can 
explore the viable thermal light DM parameter space above 
250-400 MeV. Fixed-target beam-dump experiments hold sig¬ 
nificant promise to probe lighter end of this range: future neu¬ 
trino facilities can improve sensitivity for DM below the pion 
threshold ED, while future electron beam-dumps with rel¬ 
atively small forward detectors extending the sensitivity for 
m x > ?n 7r /2 ll3Tl [721 . The fixed-target program to search 


for visibly decaying mediators, summarized in l?29l . will play 
an important role in testing the intermediate-mass scenario 
shown in Figure [2] 

Fixed-Target Missing Momentum Concept A particu¬ 
larly powerful probe of vector-portal light DM is to search for 
missing momentum in electron-nucleus fixed target collisions 
E7H23, based on the observation that most of the beam en¬ 
ergy in processes eN —> eNA', A' — > \X (and similarly for 
off-shell A') is typically carried by the invisible \X pair. Such 
experiments’ signal rate scales as oc e 2 (as opposed to e 4 for 
beam-dump searches). The authors have argued in ED that 
such an experiment can suppress backgrounds to < 1 event in 
10 16 electrons on target, provided the detector is capable of 
measuring the recoiling electron’s transverse momentum and 
vetoing products of rare photo-nuclear reactions. Figures |T] 
and [2] illustrate the significant potential of this approach. 

Scenarios Without a Thermal-Relic Milestone For DM 
heavier than the mediator, /-channel annihilation into two me¬ 
diators typically dominates over the s-channel processes used 
to compute the “relic density” curve in Figures |T] and [2] For 
scalar DM and/or vector mediators, this process is never p- 
wave suppressed, and in the inelastic scenario it allows two of 
the lighter particles to annihilate one another — thus, ther¬ 
mal DM models of these types are inconsistent with the CMB 
constraint. However, there are several scenarios compatible 
with the CMB where annihilation into two mediators dom¬ 
inates: annihilation of fermion DM into scalar mediators ip- 
wave), annihilation into vector mediators heavier than the DM 
(relying on the Boltzmann tail ll30l ). and asymmetric DM. The 
thermal relic abundance cannot be used to predict a lower- 
bound on y in these cases. A more appropriate figure of merit 
than y for comparing direct detection and accelerator-based 
experiments for mA' < m tpM * s oiane 2 (or similarly for 
scalar mediators), where the direct detection limits are eval¬ 
uated conservatively for ss m v ( x ). 

Summary This article shows that many of the simplest 
models for sub-GeV thermal DM are consistent with all cur¬ 
rent data. A generic and simple possibility is that DM couples 
to the SM through a kinetically mixed dark photon. Although 
symmetric Dirac fermion DM annihilating through the vector 
portal is excluded by CMB limits, scalar, pseudo-Dirac, and 
asymmetric DM scenarios are all largely untested. These sce¬ 
narios define sharp milestones for future experiments to reach. 
Together, planned //-factories (Belle-II, if equipped with a 
mono-photon trigger), direct detection experiments (Super- 
CDMS), and possible electron fixed-target experiments based 
on missing momentum should be capable of reaching this tar¬ 
get over almost all of the MeV-to-GeV mass range. 
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